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Abstract 

We introduce and study a magnetocaloritronic circuit element based on a domain wall that can move under applied 
voltage, magnetic field and temperature gradient. We draw analogies between the Carnot machines and possible 
devices employing such circuit element. We further point out the parallels between the operational principles of 
thermoelectric and magnetocaloritronic cooling and power generation and also introduce a magnetocaloritronic figure 
of merit. Even though the magnetocaloritronic figure of merit turns out to be very small for transition-metal based 
magnets, we speculate that larger numbers may be expected in ferromagnetic insulators. 
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1. Introduction 

There have been numerous realizations of Carnot 
machines since they were first envisaged by Nicolas 
Leonard Sadi Carnot, in both direct (i.e., engine) and 
reverse (i.e., refrigerator or heat pump) modes of op- 
eration. While the traditional mechanical Carnot ma- 
chines are based on the alternating adiabatic and isother- 
mal processes controlled by the conjugate pair of vari- 
ables (P,V), same thermodynamic principles can be put 
to work in the realizations of magnetic machines rely- 
ing on the magnetocaloric effect fljj] . These latter ma- 
chines are operating in the space of the conjugate vari- 
ables (H,M). At nanoscale, however, one might need 
to rely on different principles and thermoelectric cool- 
ing and power generation appear to be very promising 
lUHt]. In particular, large values of the thermoelectric 
figure of merit have recently been suggested for molec- 
ular junctions In this paper, we envision an inter- 
play of the magnetocaloric and thermoelectric function- 
alities. 

A growing interest in spin caloritronics that com- 
prises the spin related phenomena with thermoelectric 
effects has been spurred recently by many promising ap- 
plications ]5|, la |7[] . Thermoelectric spin transfer relates 
the heat current to magnetization dynamics |8, 9, m 
while opposite effect of heat currents resulting from 
magnetization dynamics should also occur IllOi 1 1111 . 
The spin-transfer tor q ue | l2l 1 1 311 in spin valves and 
domain walls Il5l 116] has been well understood 
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for transition-metal based magnets fn, 18, T^] which 
already led to many applications [20]. The recipro- 
cal effect to the spin-transfer torque results in electro- 
motive forces induced by the magnetization dynamics 
112 lL l22l I23L 12411 . All these pave the way for novel 



devices that can output as well as be controlled by 
temperature gradients, electric currents, and magnetic 
fields. These machines can have similar functionalities 
as Carnot machines and work at nanoscale. 



In this paper, we introduce and describe a magne- 
tocaloritronic circuit element utilizing magnetic domain 
wall motion. Further, we use this element to demon- 
strate the principle of magnetocaloritronic cooling and 
power generation. We also draw parallels between 
the operational principles of thermoelectric and magne- 
tocaloritronic cooling and power generation. This pro- 
gram naturally leads us to the introduction of the mag- 
netocaloritronic figure of merit TZmc which encodes 
information about the maximum efficiency of such de- 
vices. Our estimates indicate a very small figure of 
merit for typical transition-metal based magnets; how- 
ever, we speculate that one can achieve better efficien- 
cies using the ferromagnetic insulators in which the heat 
transferred by spin waves will better couple to the tex- 
ture dynamics in the absence of dissipation related to 
the electron-hole continuum. 
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Figure 1 : (Color online) A domain-wall based circuit element can be 
controlled by applying voltage AU, magnetic field H and temperature 
gradient AT. Here, we consider transverse head-to-head Neel domain 
wall parallel to the y axis in the easy xy plane. The constants K and 
K ± describe the easy axis and easy plane anisotropy. 



2. Magnetocaloritronic circuit element 

In order to explore new functionalities, we introduce 
a domain-wall based circuit element (see Fig. Q]) that 
combines the capabilities of a thermoelectric contact, 
heat pump and generator of electromagnetic field. The 
functionalities of this circuit element can be controlled 
directly by applying magnetic field which leads to do- 
main wall motion along the magnetic field in the direc- 
tion of the lower free energy. Alternatively, one can 
control the domain wall by applying voltage and tem- 
perature gradients which also couple to the domain wall 
motion through the viscous interaction of the charge and 
energy currents with the magnetization dynamics. The 
velocity of the domain wall then becomes: 



v = 



yHW - pfij/s - p0j q ls 



(1) 



where the domain wall acquires some velocity in re- 
sponse to the magnetic field H, the charge current j and 
the energy current j q . Equation (Q3 is derived below 
in this section from the Landau-Lifshitz-Gilbert (LLG) 
equation for the Walker ansatz along with the introduc- 
tion of the coupling constants p, p ,/3 and/? . The other 
parameters are the Gilbert damping constant a, the do- 
main wall width W, the spin density s where sm = M/y, 
with M being the magnetization density, m the unit vec- 
tor along the spin density and y the gyromagnetic ratio 
(y < for electrons). 

By writing the equation for the entropy production 
in the form analogous to the microscopic form in Ref. 
C3: 



S = 
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(2) 



and identifying the thermodynamics variables j q , j and 
X, we can phenomenologically generalize Eq. (Q]) to 
systems with the most general macroscopic equations 



describing domain wall dynamics ll 111 : 
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where X is the position of the domain wall, T = (Ti + 
T2)/2 and the last 2 equations have been inverted with 
respect to j and j q for convenience of the following 
derivations. We introduced 5 kinetic coefficients Ox, 
Oj, O/j, Oxt and Oxfi in accordance with the Onsager 
principle. 

We will now turn to the microscopic derivation of 
Eqs. (0 for the case of transition-metal based magnets 
keeping in mind that, in principle, these equations can 
be valid for other systems as well (e.g. for magnetic 
semiconductors in Ref. 12511 ). The dynamics of the cir- 
cuit element in Fig. Q] can be conveniently described 
by the following generalization to the Landau-Lifshitz- 
Gilbert (LLG) equation JT^]: 

dxP = -gj + £jq + p (m x d x m + pd x m) m , (4) 
d x T/T = %j - {jq +p'(mx d x m + 08 x m) m , (5) 

s(l + amx)m + m x H = p (d x m + fim x d x m) j 
+p (d x m +0mx <9 v m) j q , 

(6) 

where j is the charge current and j q = ju - P j is the 
energy current offset by the energy corresponding to 
the chemical potential p with ju being the ordinary en- 
ergy current. The kinetic coefficients g, f and or can 
in general also depend on temperature and texture in 
isotropic materials, for the latter, to the leading order, 
as g = go + T] g (d x m) 2 , % = £ + ^(3 x m) 2 , etc.. The 
coefficients p and p' describe the so called nondissi- 



pative 11261 1 1 OTI coupling of the magnetization dynam- 
ics to the charge and energy currents. The correspond- 
ing viscous corrections due to electron spin's mistrack- 
ing of the magnetic texture are described by the coef- 
ficients (3 and 1 3, 1^1 . The coefficients g, % and £ 
can be related to the thermal conductivity k - 1/£T, 
the Peltier coefficient II = g/g and the conductivity 
l/o" = g - T1 2 /kT. In general, H is different from the 
usual "effective field" corresponding to the variation of 
the Landau free-energy functional F [m,p, T] with re- 
spect to m at a fixed p and T, and can be expanded phe- 
nomenologically in terms of small d x T and d x p lllOll . In 
order to avoid unnecessary complications, we assume 
here that even in an out-of-equilibrium situation, when 



d x T + and d x fi + 0, H depends only on the instan- 
taneous texture mix) so that H = d m F. The texture 
corrections in Eqs. <(4j and (0 modify the energy and 
charge flows and can be relatively large in some cases 
A 1 CHI leading to texture corrections of the Gilbert damp- 
ing 12611 in the LLG Eq. (0. However, for sufficiently 
smooth domain walls, these corrections to the Gilbert 
damping are small and will be disregarded hereafter. 
The parameters p and p can be approximated in the 
strong exchange limit as ifioll : 



h o- (i - s 
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(7) 



where kq, <jq and IIo are the thermal and ordinary con- 
ductivities and the Peltier coefficient defined in the ab- 
sence of textures, respectively. The polarizations are de- 
fined as p = (crj - oj^/cro and p s = (IlJ - n^)/(nj + 
Iljj;) = en j /2n where II S 



(Til - flh/e is the spin 



crj + o-Q, and — e is the charge 



Peltier coefficient, ctq 
of particles (e > for electrons). 

We will describe the domain wall in Fig. [T]by the 
Walker ansatz valid for weak field and current biases 



123, |28|,ll 



ip(r, t) = <£(f), In tan 



9ir, f) _ x - X(t) 
2 = Wit) 



(8) 



where the position-dependent spherical angles (p and 
9 parametrize the magnetic configuration as m = 
(cosf?, sin 9 cos <f, sin 9 sin <p), Xit) parametrizes the net 
displacement of the wall along the x axis, and we as- 
sume that the driving forces (//, j and j q ) are not too 
strong so that the wall preserves its shape and only 
its width Wit) and out-of-plane tilt angle O(f) undergo 
small changes. By substituting the ansatz (0 in Eq. (0 
with the effective field given by 



H = iH + km x )x - K ± m z z + A V m, 



we obtain: 
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(9) 



W = 



where A is the stiffness constant and K and K ± describe 
the easy axis and easy plane anisotropics, respectively. 
The steady state solution of Eq. (0 below the Walker 
breakdown with O(f) = const and X — ut leads to the 
result in Eq. (Q]i. By comparing Eq. (f3]) with Eqs. 



Figure 2: (Color online) Magnetocaloritronic cooling can be realized 
by moving domain walls between two regions with different viscous 
coupling fS . The cyclic motion of the domain walls is maintained by 
the rotating clockwise magnetic field (along a horizontally elongated 
elliptical trajectory) in sync with the domain walls steadily circulating 
clockwise; this process is analogous to a pump which produces a dc 
heat flow between the cold (T c ) and hot (Ti,) junctions. The analogy 
to the thermoelectric cooling based on p- and n-type couples yj can 
be best seen when p l ft l = -p 2 /3 2 - 



(0 and (0, we can also write the following expressions 
for the kinetic coefficients: 
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where the coefficients g and f correspond to g and ( 
averaged over the wire. 



3. Magnetocaloritronic cooling 

The circuit element in Fig. Q] can in principle be 
used as a working body of a magnetic refrigerator ex- 
ploiting the magnetocaloric effect [1]. Realization of 
such a refrigerator, however, necessitates some sort of 
switch connecting and disconnecting the circuit element 
to reservoirs. 

At nanoscale, such switches can be hard to realize and 
in this work, we will explore a different route that has 
more parallels to thermoelectric cooling [2] in which the 
cooling effect appears at the junction formed by two dif- 
ferent conducting materials. In our case (see Fig. 0, 
the charge carriers are replaced by domain walls slid- 
ing along the wire due to rotating magnetic field which 
ensures cyclic motion of domain walls. As it can be ob- 
tained from Eqs. © and (0, the heat transfer originates 
from differences in the viscous coupling f3 in the upper 
and lower parts of the circuit in Fig. [2] which can be a 
result of different amounts of magnetic impurities in the 
corresponding parts. 

It is customary to describe the efficiency of ther- 
moelectric circuits by a material figure of merit Z = 
Y1 2 ct/T 2 k corresponding to the maximum achievable 
temperature difference ~ ZT 2 . In most circumstances, a 



dimensionless figure, ZT, is quoted and it corresponds 
to the relative temperature difference. Let us formulate 
an analog of figure of merit for magnetocaloritronic de- 
vice. We consider the case of relatively small temper- 
ature gradients T/, - T c <K T/,. In order to maintain 
the domain wall motion, we have to perform the work 
2MHUA per one pass and all this work is eventually 
dissipated, where H is the applied field along the mag- 
netic wires and 31 is the cross-section of the wire. From 
Eq. (0, we find that in the absence of temperature bias, 
the domain wall induces the heat current: 



;tr = 



Oxt 2MH 
~0~, ~ 



(10) 



For simplicity, we suppose that the charge current can 
not flow in the device in Fig. [2] (e.g. due to a break- 
ing point or upper and lower parts could be different 
p- and n-type semiconductors, whose junctions block 
the current flow). Supposing that the dissipated work 
is equally distributed between the reservoirs, the exter- 
nally induced heat flow from the cold/hot reservoir per 
one wire then becomes: 



cold T h ~ T c 

J „ ,, 
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(11) 



Because the domain wall cooling is proportional to the 
magnetic field and dissipative heating is proportional 
to the magnetic field squared, the maximum decrease 
in the junction temperature is obtained at an optimal 
magnetic field. This situation is similar to thermoelec- 
tric cooling in which the Peltier cooling is proportional 
to current and Joule heating is proportional to the cur- 
rent squared, leading to existence of an optimal current. 
From Eq. ( fTTT i. the maximum temperature difference 
that the domain wall motion can maintain becomes: 



(T h - T c ) 



max 



TO\t 
20 x O T ' 



We can finally introduce the figure of merit for the mag- 
netocaloritronic cooling: 
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where the last part of the equation is written for a 
domain wall described by microscopic Eqs. (0, (0 
and (0. Notice that TZmc > as it follows from 
the thermodynamic inequalities Ot > and Ox > 
which guarantee that the entropy production in Eqs. 



d2j is always positive. It is also worthwhile to note 
that the efficiency TZmc becomes infinite when X = 
2p 2 /3 2 /WLas which corresponds exactly to the lower 
bound of the thermal resistivity lioll . £ = T]^(d x m) 2 with 
j]^ = f3' 2 p' 2 las, averaged over the Walker ansatz. 

Above, we only considered dissipation in the lower 
wire in Fig. [2] Consideration of the upper wire does not 
change our results when the upper wire is mirror sym- 
metric to the lower with p l /3 l = -p 1 ji 1 (e.g. we do not 
see any principal contradictions in existence of negative 
0). In a different scenario when { = and 2 + 0, it 
is possible to minimize the effect of dissipation in the 
upper part by keeping the upper wire disconnected from 
the cold/hot junction most of the time apart from the 
moments when the domain wall moves through. 

It is interesting to see that a system in Fig. [2] made 
of p- and n-type semiconductors for upper and lower 
parts will have opposite p/3 in those parts leading to a 
device that can generate electrical power from a rotating 
magnetic field in accordance with Eqs. (0. 

At last, we recover the well known expressions for the 
maximum coefficient of performance (COP) |29jl writ- 
ten for the magnetocaloritronic cooling and heating: 



hot 
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where we maximized the above 
to H at a fixed temperature bias, 
is recovered when TZmc — > °°- 
Using Eq. (0, we can 
tocaloritronic figure of merit via 
ure of merit: 



VI + TZmc + 1 

(14) 

equations with respect 
The Carnot efficiency 

express the magne- 
the thermoelectric fig- 
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which gives a very small number Zmc ~ 1CT 7 Z for a 
typical domain wall in a Py wire at room temperature 
with L = 2W =100 nm. Much more pronounced cool- 
ing effect could be seen in MnSi below 30K for which 
we obtain Zmc ~ 10" 3 Z (MnSi parameters are taken 
to be the same with Ref. B101D . We speculate that one 
can achieve better efficiencies using the ferromagnetic 
insulators in which the heat transferred by spin waves 
will better couple to the texture dynamics in the absence 
of dissipation related to the electron-hole continuum. 
Large viscous /?-like coupling with domain wall motion 



L 




Figure 3: (Color online) Magnetocaloritronic power generation can 
be realized by maintaining cyclic motion of domain walls between 
two regions with different viscous coupling /? . The cyclic motion of 
the domain walls can be maintained by the temperature gradient when 
p fl has opposite sign for the upper and lower parts. Alternatively, 
when /Jj = 0, a small magnetic field H^ ac ^ can return the domain 
walls to the hot junction. The power (useful work) is extracted from a 
solenoid encircling one of the wires. 

have recently been predicted for dirty (Ga,Mn)As l25ll 
which can lead to larger TZmc according to Eq. (TTZl l. 
Note that the best materials available today for devices 
that operate near room temperature have a ZT of about 

1 t2 BH. 

4. Magnetocaloritronic power generation 

Thermoelectric devices find numerous applications 
as voltage generators [2]. In this section, continuing 
the analogy between the thermoelectric devices and the 
magnetocaloritronic devices, we will show how power 
and useful work can be generated magnetocaloritroni- 
cally. The device depicted in Fig [3] once again contains 
two regions with different viscous coupling . The tem- 
perature gradient propels domain walls in the lower part 
with (3 + as it follows from Eq. [3] while in the up- 
per part with [3 — the domain walls are inert to the 
temperature gradient and move due to the presence of 
a very small magnetic field H^^. As an alternative, 
one could also consider the mirror symmetric case with 
p 1 /3 l = -p 2 fi 2 anc l w i tn an additional solenoid encir- 
cling the upper wire. 

Let us calculate the ratio of the useful power to the 
losses (efficiency) for a magnetocaloritronic device. We 
again suppose that the charge current can not flow in the 
device in Fig. [3] (e.g. due to a breaking point or upper 
and lower parts could be different p- and n-type semi- 
conductors also blocking the current flow) and we con- 
sider the case of relatively small temperature gradients 
Ti, - T c <K 7ft. The dissipation due to domain wall mo- 
tion 2MHUH is again evenly distributed along the wire 
and between reservoirs as we make similar assumptions 
with the previous section. The losses then appear due to 



the finite heat conductivity j^ om = (T h -T c )/TLO T and 
due to domain wall motion. The motion of domain wall 
through the solenoid will lead to electromotive force in 
the solenoid U = 8nMJlNX/L where N is the number 
of loops in the solenoid and X is the average velocity of 
the domain wall that can be found from Eq. (0 with the 
magnetic field inside of the solenoid H = 4nIN/L. One 
can write the following expression for the useful power: 

UI = 2MJDCH , 

which, as in the previous section, suggests the existence 
of the optimal H for the maximum power outcome. We 
can finally recover the well known expression for the 
maximum thermoelectric efficiency of power genera- 
tion H29I1 written for the magnetocaloritronic power ef- 
ficiency. By maximizing the ratio of the power (work) 
UI to the losses (heat absorbed at the hot end given by 
Eq. (fTTT i) at a fixed temperature bias for the device in 
Fig. [3] we obtain: 

UI T h -T c VI + TZ mc - 1 

V = = . . (15) 

T h VI + TZ mc + TJT h 

As expected, the magnetocaloritronic figure of merit 
also contains information about the efficiency of the 
magnetocaloritronic device working as a power (useful 
work) generator. The Carnot efficiency is once again 
recovered when TZmc — > 00 ■ 

In the above estimates, we again considered dissipa- 
tion only in the lower wire in Fig. [3] Consideration 
of the upper wire does not change our results when the 
upper wire contains an extra solenoid and is mirror sym- 
metric to the lower with p l fi l = -p 2 P 2 - ^ n trie scenar i° 
when y8j = and 2 + 0, a small magnetic field #back 
returns the domain walls to the hot junction and we min- 
imize the effect of dissipation in the upper part by keep- 
ing the upper wire disconnected from the hot junction 
most of the time apart from the moments when the do- 
main wall moves through. 

5. Conclusions 

We introduced and described the magnetocaloritronic 
circuit element using recently formulated phenomeno- 
logical theory of thermoelectric spin transfer [ 10]. The 
velocity of domain wall in such a circuit element in re- 
sponse to magnetic field, charge current and energy flow 
is calculated. We also derive the most general form of 
phenomenological equations describing dynamics of a 
domain wall in response to the magnetic field, applied 
voltage and temperature biases. We conclude that such 



hybrid device combines functionalities of thermoelec- 
tric and spintronic devices. 

As an example of some of the possible functionali- 
ties of the circuit element, we propose a realization of 
magnetocaloritronic cooling and power generation. We 
further study the efficiency of the magnetocaloritronic 
cooling and power generation which leads us to the in- 
troduction of the magnetocaloritronic figure of merit by 
analogy to the thermoelectric figure of merit. Our esti- 
mates of the magnetocaloritronic figure of merit for Py 
and MnSi give very small numbers unusable for appli- 
cations. However, we speculate that one can achieve 
better efficiencies using the ferromagnetic insulators in 
which the heat transferred by spin waves will better cou- 
ple to the texture dynamics in the absence of dissipation 
related to the electron-hole continuum. 
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